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ABSTRACT 


The setting of alarm levels plays a vital role in a machinery condition monitoring and 
diagnostic system. In this research, two approaches to setting vibration alarm levels using 
vibration signals produced by fire pumps are presented in the time and frequency domains. 
In the time domain, the cross peak analysis (CPA) is proposed to extract the dominate peak 
points. The distribution of these cross peak points is found to have a lognormal distnbution 
and can be normalized to a Normal distribution in the VdB domain. The computed +26 
value in the VdB domain ts suggested for use as the alarm level. In the frequency domain, 
the 1/1 octave band analysis (OBA) 1s introduced. Three artificial fault simulations are con- 
ducted to compare the 1/1 octave band method with the broadband method. The results 
show that the 1/1 octave band method is more sensitive to the changes in VdB level than 
the broadband method. The computer programs to perform these two analyses are wnitten 


using MATLAB. Examples of the use of these programs are included in this report. 
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I. INTRODUCTION 


In recent years, the msing cost of machine maintenance has driven engineers to 
develop more economical and efficient methods to determine machine health and to plan 
accordingly the required preventive and corrective maintenance. The most popular 
technique in use today 1s a predictive maintenance program based on condition monitoring. 
In naval applications, condition monitoring 1s commonly achieved utilizing vibration 
measurement and analysis on-board surface ships and submarines [Refs. 1-7]. 

This research focuses on the use of vibration measurement to monitor machine health 
and to diagnose system problems which could lead to machine failure. In addition to 
providing accurate and understandable data on the machine’s current condition, a 
monitoring and diagnostic system must also limit the number of false alarms. Alarm 
threshold settings are vitally important in machine vibration diagnostics. Alarm thresholds 
set too high may result in premature machine failure caused by an undetected failure 
condition. Alarm thresholds set too low may result in frequent false alarms causing 
unnecessary system interruptions and repairs. False alarms also reduce operator confidence 
in the monitoring and diagnostic system. 

Because the optimum setting of alarm thresholds in vibration monitoring and 
diagnostic systems continues to be problematical, the goals of this research were: 


¢ To establish a statistical analysis method for setting vibration alarm levels by using 
time domain data. 


¢ To investigate the difference in broadband alarm levels between the cross peak 
method and vibrometer readings. 


¢ To perform a 1/1 octave band analysis for setting vibration alarm levels by using 
time and frequency domain data. 


Chapter II provides the background information on the machinery utilized in this 


research and a discussion of alarm level setting. Chapter III describes the procedure of data 


acquisition and processing and the methods of analysis. In the time domain, a statistical 
analysis method using “cross peak” data is introduced. In the frequency domain, a 1/1 
octave band concept is used to set alarm levels. The test results and observations using time 
domain data are presented in Chapter IV. Three types of shipboard fire pumps were 
analyzed using time and frequency data. Three artificial fault simulations were performed 
to assess the 1/1 octave band method. Chapter V demonstrates the 1/1 octave band analysis 
results using the NAVSSES frequency domain database. Chapter VI contains conclusions 


and recommendations for further research. 


Il. BACKGROUND 


A. DATA SOURCE 


The fire pump on a naval ship provides more than fire fighting water to damage control 
systems. The seawater provided by fire pumps is also used by vital air conditioning and 
chill water systems, main drainage systems and decontamination sprinkler systems. Atleast 
one fire pump is running whenever a ship is underway or at anchor. In this research, 23 fire 


pumps from four classes of ships are grouped into three types in Table 1. 


TABLE 1: CLASSIFICATION OF FIRE PUMP DATA 


Tyne Date of | Number of Number of Number of 
YP Survey | Fire Pumps = Locations wei aaa Sets 
10/30/89 | 6 


ee ta 
aa 
Panos [eC 







The vibration velocity signals were measured by means of transducers strategically 
placed at pickup locations on the pumps. Schematic layouts of each fire pump and pickup 
locations are depicted in Figures 1 through 3. A list of abbreviation for pickup locations 1s 
shown in Table 2. The transducer pickup placement is uniaxial with one radial pickup at 
each bearing and one axial pickup at the thrust bearing. This arrangement is listed in the 
Vibration Test and Analysis Guide (VTAG) [Ref.8]. This guide provides the most current 
technical information for each machine in a particular ship class and includes information 


on pickup locations, operating conditions and a table of vibration source components. 


Table 3 through Table 5 summanze the vibration source components. These tables identify 
exciting elements within each machine and list the vibration frequencies generated by each 
element. The vibration frequencies were normalized as multiples of the machine’s rotation 
rate, called “orders”. The time waveforms are recorded on magnetic tape by a frequency 
modulated (FM) recorder. Each pickup location on the pumps was recorded for a one 
minute time series record. Some of the important specifications of the tape recorder are 
provided in Table 6. 

In addition to the data tapes, six survey dates of data in ASCII form were received from 
the NAVSSES database. These data contain the frequency spectrum levels measured by 


vibrometer. 


TABLE 2: ABBREVIATIONS OF PICKUP LOCATION 
Abbreviation Description 
MB(FE) Motor Bearing (Free End) 
MB(CE) Motor Bearing (Coupling End) 
MB(CE/A) Motor Bearing (Coupling End/Axial direction) 
PB(FE) Pump Bearing (Free End) 
PB(CE) Pump Bearing (Coupling End) 


PB(CE/A) Pump Bearing (Coupling End/Axial direction) 
UMB Upper Motor Bearing 


LMB Lower Motor Bearing 
LMB(A) Lower Motor Bearing (Axial direction) 
UPB Upper Pump Bearing 
LPB Lower Pump Bearing 





PB(CE) PB(FE) 


is 


ate 


tH 


Pump 


PB(FE/A) 





Figure 1 Schematic Layout of a Type I Fire Pump 


TABLE 3: VIBRATION SOURCE COMPONENTS OF TYPE I FIRE PUMP 


Driver (Motor) Driven (Pump) 


Mowrsiaaty | | 1 [Poms || 
fnBindng | >| 3 [imoetervams | 6 | 6 


Operation Condition of Motor is 3595 RPM (=59.9 Hz) 
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Figure 2 Schematic Layout of a Type II Fire Pump 


TABLE 4: VIBRATION SOURCE COMPONENTS OF TYPE II FIRE PUMP 


Driver (Motor) Driven (Pump) 


owes et) |_| 1 | Pmpsan || 


Operation Condition of Motor is 3555 RPM (=59.3 Hz) 








Figure 3 Layout of a Type III Fire Pump 


TABLE 5: VIBRATION SOURCE COMPONENTS OF TYPE III FIRE PUMP 








Driver (Motor) | Driven (Pump) 


Mowrsia ety | | [Pumpman] 
KF 


Operation Condition of Motor is 3570 RPM (=59.5 Hz) 





TABLE 6: SPECIFICATIONS OF FM RECORDER 


Frequency Response 
Signal to Noise Ratio 
Harmonic Distortion 


















B. Alarm Level 

The fundamental question that must be answered before vibration condition 
monitoring can be used as a diagnostic tool for machinery is what alarm levels identify 
operation-limiting faults. Generally, the vibrations of a system can be characterized by a 
reduced data set in various domains. The cnteria used to set alarm levels can be considered 


either in the time domain or the frequency domain, as discussed below. 
1. Time Domain Criteria 


a. Vibration Severity Criterion Method 


The simplest time domain method ts the vibration seventy cnterion. The root 
mean square (RMS) value of vibration velocity is usually measured and compared with 
vibration severity charts [Ref. 9]. Figure 4 shows the vibration severity chart that was first 
introduced by T. C. Rathbone in 1939 and then later refined by the Instrument Research and 
Development Corporation. Various national standard organizations have also published 
standards for judging vibration severity. For example, Intemational Standards Organization 


(ISO) standards 2372 and 3945 provide severity guidelines for machinery. 
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Figure 4 General Machinery Vibration Severity Chart [Ref. 9] 


The disadvantages in using these criteria are (Ref. 10]: 


¢ The criteria can be used only for specific types of machinery at a standard operating 
condition. 


¢ No diagnostic information is provided. 
¢ The criteria are less sensitive during the early stages of damage. 
Thus, the vibration seventy criterion method can be only used as a rough 


indicator of machinery health. 


b. Amplitude Probability Criterion Method 
A More sophisticated method was developed using statistical analysis to 
examine the distribution of vibration amplitudes. The amplitudes used can be either peak- 
to-peak, peak or RMS readings of displacement, velocity or acceleration. Both Campel 
[Ref. 11] and Murphy [Ref. 12] use this statistical method to establish alarm levels based 
on the mean of the reading plus 3 standard deviations. The key drawback of this approach 
is that it assumes that a Normal distnbution of the linear readings exists. Vibration readings 


have a “skewed” rather than “Normal” distribution [Ref. 13]. 
2. Frequency Domain Criteria 


a. Broadband Criterion Method 


The broadband cntenion method utilizes a vibrometer which can add all the 
energy dissipated over a wide frequency range (typically 10 to 10,000 Hz). The overall 
energy is typically calculated by applying the RMS summation method to the spectrum. If 
the overall energy level exceeds a predetermined level, then alarms are tnggered. This is a 
simple and effective method in most cases. But in some cases may be inadequate because 
it is less sensitive to small changes associated with bearing defects in the presence of large 
amplitude frequency components. This limitation can be overcome by using the octave 


band or narrow band method. 


b. Octave Band Criterion Method 


The octave band cnterion method is often used in acoustics to determine the 
energy level changes due to noise and vibration. This method utilizes a constant percentage 
bandwidth to divide the frequency range of interest into several bands which provide more 
detailed information than the broadband presentation. A commonly used bandwidth is the 
one-third octave band. Early researchers used this method to check the change in each band 
level and determine if the amplitudes exceeded normal values [Ref. 14-15]. Chapter II 


discusses this method in detail. 


c. Narrowband Criterion Method 

Since the broadband and the octave band criteria methods lack detailed 
vibration information, a narrowband criterion method is gaining popularity. The bandwidth 
may be up to 10% of the center frequency range (typically 10 to 1,000 Hz). The improved 
resolution 1s generally up to 400 or 800 lines over the frequency range of interest. With 
knowledge of the forcing frequencies of rotating components, the narrowband data Is very 
useful for diagnosing specific faults. This diagnosis can be time consuming; it is also often 
difficult to explain the source of some of the peaks in the spectra. 

Frequency criteria are used in U.S. Navy surface ship to establish alarm 
settings for machinery [Ref. 3]. Vibration data for each individual machine and location are 
processed a yield mean and standard deviation values from frequency spectra data. 
Generally, broadband levels are used as a screening tools, and narrowband levels are used 


as diagnostic tools. 


1] 


II. METHODS OF ANALYSIS 


A. DATA ACQUISITION AND PROCESSING SYSTEM 


A block diagram of the data acquisition and analysis system used is depicted in Figure 
5. The Model 4800 FM recorder plays back the machinery vibration data tapes to generate 
an analog signal. By using a connector, this analog signal is distributed to the data 


acquisition software and the oscilloscope simultaneously. The oscilloscope controls the 


quality of the data by monitoring the signal time waveform. The EASYEST LX software, 
developed by Keithley Asyst, was used as a data acquisition system. This PC-based data 
acquisition system allowed us to digitize analog signals, convert the units from voltage to 
velocity and store the velocity data sequences into the hard disk of a personal computer. 
Finally, MATLAB software was used to retrieve these data sequences and perform data 
processing. MATLAB is a very powerful software tool for signal processing. All programs 
used in this research were wntten using MATLAB. The synopsis, on-line help 
documentation and program code used can be found in Appendix A. Examples for the use 
of these programs to perform the analysis are included in Appendix B. 

All of the data sets were sampled at a sampling frequency of 10 kHz. The sampling 


duration for each data set was 36.684 seconds. 


Tape Recorder 


Dallas Instrument INC. 
Model 4800 


PC-Based Data 
Acquisition Software Oscilloscope 


EASYEST-LX 


Data Processing Station 


PC 486-50MHz 


Data Processing Software 
MATLAB 


Printer 
HP LaserJet IIP 





Figure 5 Block Diagram of Data Acquisition And Processing System 


B. TIME DOMAIN - CROSS PEAK ANALYSIS 


It is obvious that the peak envelope distnbution of vibration time domain signals is 
relevant to the vibration severity of machinery component. Vibration severity criteria 
which uses measured RMS values to determine the severity of component vibration 
without considering the dispersion of the signal is truly a rough guess. A better approach is 
to represent the severity of component vibration in terms of a percentage acceptance level. 
The percentage acceptance level provides the percentage of outcomes which will not 
exceed this level threshold. 

For a Normal distribution, the acceptance level is closely related to the mean (1) and 
standard deviation (0). For example, +1.960 corresponds to a 97.5% acceptance level 
which means only a 2.5% probability the signal amplitude will exceed the u+1.960 value. 
For the sake of simplicity, the most popular assumption to assume the probability 
distribution of the signal was a Gaussian (or Normal) distnbution. However, for rotational 
machinery vibration, since we are only concern with the magnitude of the signal, the 
probability density function has a skewed rather than a Normal distribution. At this point, 


itis worthwhile to review the statistical moments and central moments of a random variable 


{Ref.16]. For the discrete random variable x with probability Pr(éx<), the nth moment E(x”) 


is defined as: 


E(x") = ya Perce) (1) 


k=] 


The first moment is particularly important in many applications and is given the name mean 


value (1): 


u, = E(x) = yx, PG (2) 


k=] 


Of greater significance are the central moments which defined as: 


a 


El (x-u)") = 3 -H)” Prep (3) 


k=] 


The central moment for 7=1 is zero. The central moment for 77=2 is a very important which 
quantity called the vanance of the random variable; 


00 


El(x-w)7) = YG - 4)? Prep (4) 


k=] 
The standard deviation, which corresponds to the dispersion of the random vanable, can be 


obtained by taking the square root of the vanance,i.e., 


o = Variance (5) 


As noted earlier, the mean and standard deviation are often used to compute the acceptance 
level. 

The higher order central moments (n>2) are often normalized by dividing by the 7th 
power of the standard deviation. The third and fourth normalized central moments are 
primarily used to indicate the shape of the probability density function. They are defined 


by the following equations: 


E[(x-p,)° 
Skewness = = Salad (6) 
o 


Kurtosis = ae (7) 
GO 


The skewness shows information about the position of the peak density relative to the mean 
value. The kurtosis indicate the spread in the distribution. For a perfect Normal distribution, 
the skewness is zero and the kurtosis 1s three. 


Since we are interested in the peak envelope distnbution of the vibration signal, we 
use the cross peak data points instead of the overall peak data points. The cross peak 


technique can extract the dominant component from a complicated time waveform data set 


by selecting the maximum peak value between every zero crossing. Figure 6 illustrates the 
difference between overall peak and cross peak points. 

As we will see in Chapter IV, the probability distributions of cross peak points have 
highly positive skewness and can be approximated as a lognormal distribution. This means 
that the distributions will become normal if the data 1s transformed by a log or natural log 
algonthm. The method chosen to accomplish this transformation was to convert the linear 
velocity reading to velocity dB (VdB) reading. The VdB 1s defined as: 

Velocity(VdB) = 20 log — (8) 
ref 
where V,.is the reference level, normally OVdB=10°° m/sec. 


Figure 7 shows a flow chart for our cross peak analysis technique. First, the sampled 
data is imported into a MATLAB worksheet. The DC offset is then removed by subtracting 
the mean value of the sampled data. The cross peak points between every zero crossing is 
then found using a subroutine. As can be seen from Table 1, the dynamic range of the tape 
recorder is 40 dB. This means that amplitude smaller than 1% of the maximum peak may 
be contaminated by noise. For this reason, we set a 1% threshold to eliminate these 


contaminated points. Equation 8 was then used to transform the linear velocity scale into 
the VdB domain. The reference 0VdB was 10°8 m/sec. After we computed the mean, 
standard deviation, skewness and kurtosis, plots of the probability density function were 


generated. Finally, the statistical analysis report can be sent to laser printer. 
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Figure 6 Comparison of Overall Peak and Cross Peak Points (a) Overall Peak Points 


(b) Cross Peak Points 
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Figure 7 Cross Peak Analysis Flow Chart 


C. FREQUENCY DOMAIN - 1/1 OCTAVE BAND ANALYSIS 
This method divides the frequency spectrum into constant percentage bands having 
the same ratio of bandwidth to center frequency [Ref. 17]. Each band has an upper 


frequency limit (/>) and lower frequency limit (77). The center frequency (/,.) of any band is 


defined as 


the = J l hy (Y) 
The ratio of center frequencies of successive proportional bands 1s the same as />//; for 


any one band. 1.e., 


= 
Je = J, = ¥2 = 9”? (10) 
fi fe WS 
One third of an octave is defined as n=1/3. The American National Standards Institute 
(ANSI) preferred center frequencies and pass bands for 1/3 octave bands ts tabulated at 
Table 7. 
The method chosen to examine the frequency domain alarm levels in this research was 
to use 1/1 octave bands (#=1). For the 10-5,000 Hz frequency range of our system, the 
frequency spectrum can be divided into 9 bands using the ANSI preferred center 


frequencies. The center frequencies and pass bands covering the frequency range 10-5,000 


Hz in 1/1 octave bands are given in Table 8. 


TABLE 7: ANSI PREFERRED CENTER FREQUENCIES AND PASS BANDS 
FOR 1/3 OCTAVE BAND 
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TABLE 8: ANSI PREFERRED CENTER FREQUENCIES AND PASS BANDS 
FOR 1/1 OCTAVE BAND 


Center Pass Band (Hz) 
Frequency (#2) Upper Limit 










Figure 8 shows the flow chart for our 1/1 octave band analysis using the MATLAB 
program developed. This technique used frequency domain data to compute the 1/1 octave 
band levels. In order to have better resolution at lower frequency, it was necessary to use 
two spectrum data with different sampling frequencies. For the time domain data, the lower 
sampling frequency spanned 2 kHz and the higher sampling frequency spanned 10 kHz. 
Twenty blocks of data is sampled at each sampling frequency, a block of data being 2048 
data points. After sampling, a FFT was performed for each block with a Hanning window 
to obtain the smoothed linear velocity spectrum. The averaging VdB spectrum was then 
produced by transforming the linear spectrum to the VdB domain and taking an average 
over the 20 data blocks. For the frequency domain data received from the NAVSSES 
database, 10 order and 100 order data sets were used as the lower frequency spectrum and 
higher frequency spectrum, respectively. The Ist through 6th 1/1 octave band levels were 
then computed using the lower frequency sampling spectrum and the 7th through 9th 1/1 
octave band levels obtained using the higher frequency sampling spectrum. Finally, the 1/ 


1 octave band analysis was done after combining these two ranges into 9 band levels. 
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Figure 8 Flow Chart of 1/1 Octave Band Analysis Technique 


IV. RESULTS OF USING TIME WAVEFORM TAPES 


A. CROSS PEAK ANALYSIS 


1. Probability Distribution of Cross Peaks 

The available data was digitized and then grouped into three types after 
examination of the spectrum patterns. Typical time sens and frequency spectrums for each 
type of fire pump are depicted in Figure 9 through Figure 11. The data collection sheets are 
summarized in Table 9 through Table 13 for each class of ship. One set of data is missing; 
the actual number of data sets analyzed was 134. The operation speed of the shaft was 
selected as a reference speed. This speed was used to normalize the frequency scale into 
“orders”. The “B.B. Level” is the broadband level measured by vibrometer when the data 
was recorded. This level represents the total energy of the velocity reading for each 
measurement. The full scale represent the calibration factor. The “B.B. Alm Level” is the 
broadband alarm level set by NAVSSES using the vibrometer reading data. Both the “B.B. 
Level” and “B.B. Alm Level” was measured for a frequency range from 10 Hz to 10,000 
Hz. The remarks column describes the quality of the data set. 

In order to understand the distributions of the cross peak points, probability 
histograms for each set of data were plotted using linear velocity scale. Figure 12 and 
Figure 13 give an example of probability histograms for type I fire pumps measured at 
MB(FE). Inspection of these histograms shows the distributions to have an exponential 
shape for the most part. With a linear velocity scale, this type of distnbution means that the 
population is mainly dominated by small peaks which ts an expected result for machinery 
in good condition. Those distributions having a large population at higher amplitudes are 
regarded as representing damaged pumps. Examination of these damaged data sets shows 


their broadband levels to be higher than the others and closer to the broadband alarm levels. 
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Figure 9 Typical Time Series and Frequency Spectrum of Type I Fire Pumps 
(10/30/89) Measured at MB(FE) 
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Figure 10 Typical Time Series and Frequency Spectrum of Type IJ Fire Pumps 
Measured at MB(FE) 
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Figure 11 Typical Time Series and Frequency Spectrum of Type [II Fire Pumps 
(3/10/89) Measured at UMB 
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TABLE 9: DATA COLLECTION SHEET FOR TYPE I (10/30/89) FIRE PUMPS 
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TABLE 10: DATA COLLECTION SHEET FOR TYPE I (2/26/90) FIRE PUMPS 


Machine Speed Pickup yak 
Name (RPM) Location (VaB) 


Remarks 








MBE) [38] 10 | _107 | _—_ 
BCE) | 97 | 100 | 109 [ 
a rPBCCE) [108 [130 [_is_[ 
| 

PBrEAY[ tt | 130 | ar | 

MBrE) [97 | 0_| orf 

vos [MBCEY | 96 | 120 [ 108 
#2 PBCE) | wos [0 ef 
0 | _9_ | a 

PBEAY| 108 | 130 |__| 


PMBE) [99 [130 [_1o7_-[ 
Fire Pump | ygyq [MBCCEY | _100_[ 130 [109 | 
"8 
PREY | 11s | _120_| _119__| Damaged 





rBREAS| rT | a0 a 
ree) [8 | 130 [07 | 

Fire Pump | see [MBCEY | 97] 20 [109 | 
4 PBKCE) | 1s [ 130 [ ef 
Pere) || a0 ff 

PREY] 12 | 0 | aif 

Bre) | 99 | ao | 7 

Fire Pump | yee [MBCCEY | ior_[ 130 [109 | 
ts rPBce) [109 [os 
Pare) | no [of 

PBrEAY] 3s [0 | er | 

Bre, | 99 [of orf 

Fire Pump | ygyg [MBCEY | _98 [120 [ 109 | 
#6 
Pare) [2 | of ef 

parE] 13 _[ 30 [ ar 


28 


TABLE 11: DATA COLLECTION SHEET FOR TYPE II FIRE PUMPS 


Machine Speed Pickup 
Name (RPM) Location 
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TABLE 12: DATA COLLECTION SHEET FOR TYPE ITI (11/24/89) FIRE PUMPS 
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TABLE 13: DATA COLLECTION SHEET FOR TYPE ITI (3/10/89) FIRE PUMPS 
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Probability Density Function for Cross Peak Points 
Type | Fire Pump (10/30/89) 
Pickup Location : MB(FE) 
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Figure 12 The Probability Distributions of Cross Peak Data for Type I (10/30/89) Fire 
Pumps Measured at MB(FE) 
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Probability Density Function for Cross Peak Points 
Type | Fire Pump (2/26/90) 
Pickup Location : MB(FE) 


‘ (a) Fire Pump #1 


£ 
4 
6 0.5 
2 
oO. 
0 
0 1 2 3 4 
Velocity (m/sec) x 10° 
(c) Fire Pump #3 
1 
2 
3 
= 0.5 
= 
Oo. 
0 
0 1 2 3 4 
Velocity (m/sec) x 10° 
(e) Fire Pump #5 
1 
2 
3 
0.5 
2 
oO. 
0 
0 1 2 3 4 
Velocity (m/sec) eae 


(b) Fire Pump #2 


Probability 





0 1 2 3 4 


(d) Fire Pump #4 


2 
3 
805 
© 
oO 
0 
0 1 2 3 4 
Velocity (m/sec) x 10° 
(f) Fire Pump #6 
1 
2 
3 
$05 
e 
oO. 
0 
0 1 2 3 4 
Velocity (m/sec) x 10° 


Figure 13 The Probability Distributions of Cross Peak Data for Type I (2/26/90) Fire 
Pumps Measured at MB(FE) 


2. Statistical Analysis Results 
As mentioned before, the distributions of the cross peak points in a linear velocity 
scale have an exponential shape histogram. This shape histogram can be approximated as 
a lognormal distribution. Because a velocity decibel is a unit of measurement equal to the 


logarithm of the ratio V/V_,,, it is convenient to use the VdB transformation to convert an 


exponential shaped distribution into a Normal distribution. Before we applied this 
transform, the 1% threshold was used to remove distorted data. The data sets with same 
type of fire pump and same pickup location were then added together to form a combined 
data set to perform the statistical analysis. The results are collected in Appendix C. Figure 
C.1 through Figure C.15 present the statistical analysis results for all of the available data. 
Figure C.16 through Figure C.21 show the statistical analysis results without including the 
damaged data sets. Table 14 tabulates these results in the linear velocity and VdB domains. 
A comparison of alarm thresholds between the broadband alarm levels (B.B. Alm. Level) 
and the computed |1+2o0 level is shown in Table 15. It should be mentioned that these two 
alarm thresholds are based on different methods using data from different domains. The 
broadband alarm level, as noted earlier, is obtained by using vibrometer readings in the 
frequency domain. The computed +o levels represent the 97.5% acceptance level of the 
cross peak envelopes in the VdB domain. Several observations can be made: 


¢ The distribution of the cross peak points in the linear velocity domain have a highly 
positive skew. This can be treated as a lognormal distnbution. 


¢ By using the Vdb transformation, the distributions in the linear velocity domain can 
be transformed to Normal or near Normal distributions. Generally speaking, the data sets 
at pickup locations for the motor have a skewness less than 0.1 and the data sets at pickup 
locations for the pump ends have a skewness varying between 0.1 and 0.5. 


¢ As can be seen from Table 14, removing the damaged data sets does not cause the 
probability distributions in the Vdb domain to become much closer to Normal as measured 
by the skewness and kurtosis. For the MB(CE) pickup location of a type II pump, not 
considering the damaged pump data sets makes the skewness even more severe than when 
the damaged pump data sets are considered. Since there were only 3 type II pumps, this 
result is not surprising. 


¢ In general, the B.B. Alarm Levels and computed t1+20 levels in the VdB domain are 
quite close. The relative errors between these two methods for Type I fire pumps are less 
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than 5%. It is also noted from Table 15 that the computed 1+2o levels in cases with small 
skewness is closer to the broadband alarm levels than in cases with high skewness. This is 
because the computed +20 levels for the Normal distribution is very sensitive to 
skewness. 
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TABLE 14: STATISTICAL ANALYSIS RESULTS FOR CROSS PEAK POINTS 
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B. 1/1 OCTAVE BAND ANALYSIS 


1. 1/1 Octave Band Analysis Results 


The 1/1 octave band analysis was performed for all types of fire pumps with the 


time domain data. Band 1 through band 6 were computed by using f,=2 kHz to provide 
better resolution at lower frequencies. Band 7 to band 9 were computed by using f,=10 kHz. 


Figure 14 through Figure 16 show the summanzed 1/1 octave band levels for each type of 
fire pump without considering the damaged data sets. The 10-5,000 Hz broadband levels 
are also shown on these plots. It was computed by adding all the VdB spectrum using a 
RMS algorithm. The means (1) and standard deviations (6) of each octave band levels for 
the same type fire pumps with the same pickup locations was also computed. The thin line 
represents the mean (1) band level and the thick line represents the mean plus one standard 
deviation (uU+o) band level. Table 16 tabulates these band levels. 

Generally speaking, the dominate levels are located in the first six bands (10 to 
710 Hz). This implies that the energy of vibration 1s concentrated at lower frequencies. For 
fire pumps with a 3555 RPM operation speed, these bands are approximated up to 12 
orders. Upon examination of the third band through the sixth band for each data set, it 1s 
obvious that the band levels for damaged data sets are higher than the for undamaged pump 


data sets. 
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1/1 Octave Band Analysis 
Type | Fire Pump 
Each Data Set Use 20 Averaged FFT Data 


(a) Pickup Location : MB(FE) 
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Figure 14 The Summarized 1/1 Octave Band Levels for Type I Fire Pumps 
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(c) Pickup Location : PB(CE) 
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Figure 14 The Summarized 1/1 Octave Band Levels for Type I Fire Pumps (Cont.) 
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1/1 Octave Band Analysis 
Type Il Fire Pump 
Each Data Set Use 20 Averaged FFT Data 


(a) Pickup Location : MB(FE) 
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Figure 15 The Summarized 1/1 Octave Band Levels for Type II Fire Pumps 
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(c) Pickup Location : MB(CE/A) 


Use 3 data sets 
10-5k Hz B.B. +6 Level=114.67 dB 
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Figure 15 The Summarized 1/1 Octave Band Levels for Type II Fire Pumps (Cont.) 
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1/1 Octave Band Analysis 
Type Ill Fire Pump 
Each Data Set Use 20 Averaged FFT Data 


(a) Pickup Location : UMB 
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10-5k Hz B.B. +0 Level=100.67 dB 
yw: 


VdB (OVdB=1¢-8 m/s) 





Frequency (Hz) 
(b) Pickup Location : LMB 
Use 6 data sets 


10-5k Hz B.B. pt+o Level=102.45 dB 
[he eee 


VdB (OVdB=16-8 m/s) 





Frequency (Hz) 


Figure 16 The Summarized 1/1 Octave Band Levels for Type III Fire Pumps 
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(c) Pickup Location : LMB(A) 


Use 8 data sets 
10-5k Hz B.B. p+o0 Level=101.73 dB 
i ae 


VdB (OVdB=1¢-8 m/s) 





Frequency (Hz) 


(d) Pickup Location : UPB 


VdB (OVdB=1¢e-8 m/s) 





Frequency (Hz) 


(6) Pickup Location : LPB 


Use 8 data sets 
10-5k Hz B.B. p+o0 Level=111.73 dB 
Lie 


VdB (OVdB=10-8 m/s) 





Frequency (Hz) 


Figure 16 The Summarized 1/1 Octave Band Levels for Type III Fire Pumps (Cont.) 


2 |< Is 
On 
SS oh gio 
on = WO oO a ee pa + Ww ™ 
= 06 oD S = = a 
6 = Ol<¢ 
00 = < 
“ oO eS 
D a 
00 = 
“NN 2 = pa = 
~ o ~ 
ai OO S 
s: 
o o 
00 | > 0 x an |r 
a or on on 
m~ |o 
e = 


69.4 


OVdB=le-8 m/sec 


peas es mo ms] ea] as 


Band Level (VdB) 


OCTAVE BAND LEVELS 


es 

3 
AZ 
= 
Qu 


Location 





MB(FE) 
MB(CE) 
PB(CE) 
PB(FE) 
PB(FE/A) 
MB(FE) 
MB(CE) 
MB(CE/A) 
PB(CE) 
PB(FE) 
LMB(A) 


TABLE 16: SUMMARIZED MEAN (1) AND STANDARD DEVIATION (06) OF 1/1 


2. Artificial Fault Simulations 


a. Fault Simulation 
In order to assess the 1/1 octave band method, three cases of artificial fault 
simulation were performed. Before we simulated these cases, we had to define a “fault”. 
For simplicity, we assumed a fault could be approximated with 5 bars, as shown in Figure 


17. The bandwidth of the fault is about 4.5 Hz. 


Frequency 





Figure 17 Fault with 6VdB Gain at Center Bar 


Based on Eq. (8), if the vibration amplitude doubles with a linear velocity 
scale then the VdB level will increase 6VdB at the corresponding frequency. Therefore, we 
defined a fault as having a 6VdB gain at the center bar and use four bars with 3 VGB (half 
the power of the center bar) and 2VdB at the side bands to simulate leakage effects around 


the corresponding frequency. 
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b. Simulation of A Misalignment Fault 
Misalignment occurs when the center lines of two shafts are offset or meet at 
an angle. The characteristics of misalignment in a spectrum include: 
¢ High amplitude axial peaks and radial peaks at 1, 2 and 3 orders of shaft RPM. 


¢ Higher harmonics of the shaft RPM (greater than 4 orders of shaft RPM) are 
generally low in amplitude. 


Thus, it is reasonable to add the 6VdB gain fault at the first and second order. 
Figure 18(a) and 18(b) show the spectrums before and after imposing the fault. Figure 1 8(c) 
compares the differences in the 1/1 octave band levels and the 10-5,000 Hz broad band 
level for these two conditions. As can be seen in Figure 18(c), the 1/1 octave band levels 
have a 4.46 VdB gain at third band (which corresponds to 1 order of shaft RPM) and a 3.13 
VdB gain at the forth band (which corresponds to 2 orders of shaft RPM). However, the 10- 
5,000 Hz broadband level only increased 2.17 VdB. Obviously, the 1/1 octave band method 


is more sensitive than broadband method. 
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Case |:Misalignment Fault 


Type | Fire Pump 
Pickup Location : PB(CE) 





(a) Good Condition (6) Damaged Condition 
120 120 
110 110 
5 5 
© 100 © 100 
90 90 
80 ‘= 
0 2S 10 0 5 10 
Orders Orders 


(c) 1/1 Octave Band Levels 


10-5kHz B.B. Level 


—— Good 107.96db 
Bad 109.79db 


1e-8nVs) 
nN 
© 
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Co 
O 

nN 


Frequency (Hz) 


Figure 18 Simulation Result for Artificial Misalignment At Shaft 
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c. Simulation of A Looseness Fault At Impeller 

Impeller looseness is a rotating element looseness. The important 
characteristics of looseness in a spectrum include: 

¢ Presence of a large number of harmonics. 
¢ Presence of half-harmonics. 

For example, the fire pump of a Type I has six impellers contnbuting to a 
forcing frequency of 6 orders of the shaft RPM, and we expect higher levels at 6 orders and 
its harmonics (12 orders, 24 orders,... etc.) for a loose impeller. The levels at half- 
harmonics (e.g., 3 orders and 9 orders) will also increase. Thus, in this looseness fault 
simulation, an artificial fault was applied by adding a fault with 6VdB at center bar at the 
6th order(i.e., 6 x shaft RPM) and adding two fault with 3 VdB at center bar at the 3rd order 
and the 9th order. Figure 19(a) and 19(b) show the spectrum of good and damaged fire 
pumps. Figure 19(c) compares the differences in 1/1 octave band levels and the 10-5,000 
Hz broadband level for these two conditions. At the fifth of the 1/1 octave band levels, there 
is a 1.5 VdB gain. However, the 10-5,000 Hz broadband level only increased 0.29 VdB. 


Again, the 1/1 octave band method is more sensitive than broadband method. 
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Case II:Looseness Fault At Impeller 


Type | Fire Pump 
Pickup Location : PB(FE) 








(a) Good Condition (6) Damaged Condition 
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80 
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(c) 1/1 Octave Band Levels 
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Figure 19 Simulation Result for Artificial Looseness At Impeller 
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d. Simulation of A Bearing Fault 
For a steady state condition, some periodic signatures exist which relate to 
corresponding bearing faults, called bearing frequencies. These bearing frequencies can be 
found in the VTAG. In this section, an artificial fault imposed at the motor bearing (free 
end) of a Type I fire pump has been simulated. The bearing frequencies of this bearing are 


tabulated in Table 17. 


TABLE 17: BEARING FREQUENCIES OF MRC310 BEARING 


Bearing Frequency Symbol 
Train order of rolling f 0.38 
element ; | 
Relative rotation order of £ 0.619 
rotating raceway : | 
Spin order of rolling f 1.98 
elements : | 
Irregularity order of f 3.96 
rolling element ie | 
Irregularity order of f. 4.95 
rotating raceway a | 
Irregularity order of f 3.04 
stationary raceway a | 


Suppose, for example, that there is a wear degradation in the inner raceway. 













Then the presence of the fundamental f;, tone with harmonics would be expected. Figure 


20(a) shows the spectrum of a good bearing. The 6 VdB fault is imposed at the 4.95 order 
as shown in Figure 20(b). Figure 20(c) shows the resulting increase in the 1/1 octave band 
levels. The 4.95 order is located at the fifth of the 1/1 octave bands. The 10-5,000 Hz 
broadband levels before and after damage are almost unchanged. However, the fifth octave 
band level has a 1.57 VdB increase. This illustrates why a single broadband level cannot be 


used to determine the condition of a machine. 
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Case Ill:Bearing Fault 


Type | Fire Pump (11/30/89) 
Pickup Location : MB(FE) 


(a) Good Condition (b) Damaged Condition 
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Figure 20 Simulation Result for Artificial Bearing Fault At Motor 
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V. RESULTS OF USING NAVSSES FREQUENCY DOMAIN 
DATABASE 


Since digital vibration survey instruments are more popular and easy to use than other 
types, only frequency spectrum reading were recorded and stored into the database. Thus, 
our 1/1 octave band analysis was performed using frequency domain data. Data measured 
for the Type I fire pumps in six survey dates are available in ASCII format. The total 
number of fire pumps analyzed was twenty-seven. Table 18 shows the available frequency 
data and corresponding survey date. Each set of data contains a 10 order and a 100 order 
spectra. The resolution for each set of data is 400 lines. The averaged spectra using these 
data are shown in Figure 21 and Figure 22 for each order range. The thin lined curve 
represents the mean levels and the thick lined curve represents for the mean plus one 
standard deviation levels. The 10 order spectrum was used to compute the first six 1/1 
octave band levels and the 100 orders spectrum was used to compute the seventh and above 
band levels. A summary of 1/1 octave band analysis results for each pickup location using 
all the fire pumps is shown in Figure 23. As can be seen from Figure 23, the higher 
frequency levels (seventh band to ninth band) are about 95 VdB for the motor bearing and 
110 VdB for the pump bearing. However, these higher frequency levels in Figure 14 using 
time waveform tape data are about 80 VdB for the motor bearing and 90 VdB for the pump 
bearing. As we reviewed the narrow band plot for the NAVSSES frequency domain data, 
we found that the levels below 54 VdB have been cut off and forced to 54 VdB. This is why 


the higher frequency levels using the NAVSSES frequency domain data are higher than the 
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levels from using time waveform tape data. Figure 24 shows a good example of a cut off 


spectrum for a motor bearing using the NAVSSES frequency domain data. 


Table 18: DATA COLLECTION SHEET FOR AVAILABLE FREQUENCY 
DOMAIN DATA FROM NAVSSES DATABASE 


Number of Fire | Numberof Pickup | Number of Data 
Survey Daie 
Pumps Location Set 
aun | | s e 


a 











an 











07-JUN-90 5 





an 





Narrow Band Average for 0-10 Order 
Using NAVSSES Database 
Type | Fire Pump 


Each Position Use 27 Pumps for Averaging 


(a) Pickup Location : MB(FE) 


VdB (OVdB=1e-8m/sec) 





Frequency (Orders) 


(b) Pickup Location : MB(CE) 


VdB (OVdB=1e-8m/sec) 





Frequency (Orders) 


Figure 21 The Averaged Spectra for Type I Fire Pumps (0-10 Order) Using 
Frequency Domain Data 
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(c) Pickup Location : PB(CE) 
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(d) Pickup Location : PB(FE) 
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(d) Pickup Location : PB(FE/A) 
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Figure 21 The Averaged Spectra for Type I Fire Pumps (0-10 Order) Using 
Frequency Domain Data (Cont.) 
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Narrow Band Average for 0-100 Order 
Using NAVSSES Database 
Type I Fire Pump 


Each Position Use 27 Pumps for Averaging 


(a) Pickup Location : MB(FE) 
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(b) Pickup Location : MB(CE) 
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Figure 22 The Averaged Spectra for Type I Fire Pumps (0-100 Order) Using 
Frequency Domain Data 
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(c) Pickup Location : PB(CE) 
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Figure 22 The Averaged Spectra for Type I Fire Pumps (0-100 Order) Using 
Frequency Domain Data (Cont.) 
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1/1 Octave Band Analysis 
Using NAVSSES Database 
Type I Fire Pump 


Each Position Use 27 Pumps 


(a) Pickup Location : MB(FE) 


100 order B.B. u+o0 Level=106.85 dB 
eee 


VdB (OVdB=16-8 m/s) 
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(b) Pickup Location : MB(CE) 


100 order B.B. u+o Level=111.07 dB 
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Figure 23 The Summarized 1/1 Octave Band Levels for Type I Fire Pumps Using 
NAVSSES Frequency Domain Data 
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(c) Pickup Location : PB(CE) 
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Figure 23 The Summarized 1/1 Octave Band Levels for Type I Fire Pumps Using 
NAVSSES Frequency Domain Data (Cont.) 
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(d) 0-100 order for MB(CE) 
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Figure 24 The Narrow Band Spectra for Type I Fire Pumps Using NAVSSES 
Frequency Domain Data Measured at MB(FE) and MB(CE) 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


Both time and frequency domain analyses were performed using fire pump vibration 
data to determine an appropriate alert level. Cross Peak Analysis in the time domain and 1/ 
1 Octave Band Analysis in the frequency domain were used to set the alert threshold 
level(s). The following conclusions are drawn: 

Time Domain Analysis (Cross Peak Analysis) 


a. The measured peak envelop data at five different pickup locations follows a 
Gaussian probability distnbution in the VdB domain well. 


b. The +20 value computed using a peak envelop probability density function in the 
VdB domain gives a broadband peak amplitude alert level rather than the energy content in 
the vibration signal. 


c. We used 12 measured data sets (6 fire pumps and 2 different measurement dates) 
which were available on cassette data tapes. To improve the quality of averaging, more 
measured data sets would be required. 


d. NAVSSES does not save time domain data on any storage media at this time. 
However, for some shipboard machinery, such as low and high pressure air compressors, 
time domain data is required to cast measured data in the time-frequency domain to detect 
the possible faults. 
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Frequency Domain Analysis (1/1 Octave Band Analysis) 


a. 1/1 Octave Band Analysis (OBA) uses nine frequency band alert levels which 
provides more detailed information about machine condition than the simple broadband 
level measured by a vibrometer. When a significant change in VdB occurs in a particular 
frequency band or bands, narrow band zoom mode analysis can be performed for the 
selected frequency band(s) to identify the component(s) which may have faults. 


b. OBA divides the frequency range into 10 bins over 10 kHz. The VdB level in each 
frequency bin is quite sensitive to changes in the energy content of the measured vibration 
signal. 


c. The OBA approach can be easily incorporated into the Digital Vibration Survey 
Instrument (DVSIJ). 


d. Five sets (27 fire pumps) of 400 line frequency domain DVS] data were received 
from NAVSSES and OBA was performed. The VdB levels in the higher frequency bins 
were generally high. This was caused by the fact that the vibration signals which were 
lower than the VdB limit(54 VdB) were set to the lower VdB limit in a specified dynamic 
range. 


e. OBA was also performed using time domain data obtained from cassette data tapes. 
The changes in VdB levels at a particular frequency and the corresponding sidebands were 
introduced in narrow band FFT. The results are transferred into the octave band domain and 
the sensitivity in VdB levels in frequency bins was analyzed. The results are quite 
promising. 


f. NAVSSES maintains 400 line frequency domain data for all machineries in a 


graphics postscript format, not in ASCII or digitized format. It was not possible to convert 
this graphics format data to ASCII format. 
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Based on current studies, the following recommendations are made: 


a. In the current machinery monitoring practice, frequency domain data is commonly 
saved using DVSI. However, for some shipboard machinery such as low and high pressure 
air compressors, time domain data is required to cast the measured data into the time- 
frequency domain to detect the possible faults. We suggest saving time domain data for the 
limited number of shipboard machinery. 


b. The 1/1 Octave Band Analysis Method is suggested to implement in DVSI. 


c. NAVSSES maintains 400 line frequency domain data for all machinery in graphics 
postscript format, not in ASCII or digitized format. It was not possible to convert this 
graphics format data to ASCII format. We suggest saving future files in ASCH format to 
permit further analysis of the data. 


d. It is suggested to saving/storing DVSI data “as is”, without alteration. 


APPENDIX A. MATLAB PROGRAM CODE 


A. ON-LINE HELP DOCUMENTATION 


This section of the guide contains all of the MATLAB functions that we have 
developed during this research. Each function contains a purpose, synopsis and description 
in alphabetical order. The MATLAB built-in functions used in these programs are 
described in the MATLAB reference guide. 

The code was written for the MATLAB 4.0 version. For those who want to use these 
programs in older versions need to consult the MATLAB reference guide for necessary 


modifications. 
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1. kurt 


Purpose 
kurtosis or fourth moment of vectors and matrices. 
Synopsis 
kurtosis=kurt(x) 
D — 
krut calculates the kurtosis or fouth moment value. 
For vectors, kurt(x) is the kurtosis value of the elements in vector x. For 
matrices, kurt(x) is a row vector containing the kurtosis value of each 


column. 
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2. logbar 
Purpose 


Plot 1/1 octave band bar plot 

Synopsis 
logbar(Vdb) 
logbar(Vdb,LineStyle) 
logbar(Vdb,LineStyle,LineWidth) 

D =i 
logbar(Vdb) draws a 1/1 octave band bar plot using Vdb as the y axis 
and the ANSI preferred center frequencies and pass bands as the x axis. 
LineStyle and LineWidth are optinal parameters. The default LineStyle 
is a solid line. See “plot” for more information about LineStyle. The 


default LineWidth is 0.5. 
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3. mstat 


Purpose 


Statistical analysis for time domain data. 


Synopsis 


mstat(y,fs,nbar,id) 


mstat(y,fs,nbar,id,manaly,mid,locat,mregen,pn) 


pene 


mstat(y,fs,nbar,mid,locat,manaly,mthre,z,mregen,pn) will analyze the 


original data y using the manaly method and produce a statistical 


report with a nbar pdf plot. 


mstat(y) allows the user to execute the program in an interactive way. 


y 
fs 


nbar 
mid 
locat 


manaly 


mthre 


z 


mregen 


pri 


amplitude (calibrated velocity data) 
sampling frequency (Hz) 

number of bars in Probability. Density Distnbution Plot 
type of machinery 

pickup location 

method of analysis 

=] means All data points analysis 

=2 means Overall peak points analysis 
=3 means Cross peak points analysis 
with threshold or not 

=1 without threshold 

=2 with threshold 

percentage of threshold 

transform to Vdb or not 

=] without transform 

=2 transform to Vdb domain 


print or not (‘y’es or ‘n’o) 
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npc 
Purpose 
Cross peak points. 
Synopsis 
ype=npc(y) 
D — 
ycp=npc(y) will find the peak points between each two zero-crossing.. 
y original data 


ycp cross peak data 
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5. npeak 
Purpose 


Overall peak points 
Symopsis 
yp=npeak(y) 
eye 
yp=mpeak(y) will find the overall peak points. 
y original data 


yp overall peak data 
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6. 


nthre 


Purpose 


Threshold 


Synopsis 


yt=nthre(y,z) 


Peccenis 


yt=nthre(y,z) will remove the points with z% threshold . 


a 
yt 


Z 


original data 

data after threshold 

percentage of the threshold 

If z=1 then the program will remove the peak data small than 


1% of the maxinum peak value. 
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ce 


octave 


Purpose 
1/1 octave band level. 
Synopsis 
oct=octave(y) 
oct=octave(y,1),12) 
D wae 
oct=OCTAVE(y) returns the 1/1 octave band Vdb values from band 1 
to band 10. 
oct=OCTAVE(y,1],12) returns the 1/1 octave band Vdb values from 


band i) to band 1, and the others bands will be assigned a value of 0. 


y(:,1) : Vdb values 
y(:,2) : frequencies 
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psd 
Purpose 


To evaluate the averaging VdB values using Hanning window. 
Synopsis 

VdB=psd(y,nblock) 
D oe 

VdB=psd(y,nblock) returns the average VdB values using time domain 


data y. The number of averagings is specified by mnblock. The 


recommend number of averagings is 20. The data matrix y must be a 2” 
by 2 matrix. y(:, /) contains time series and y(:, 2) contains corresponding 


velocity data. The hanning window is used to damp out the effects of the 


Gibbs phenomenon. The reference 0VdB is 10°8 m/sec. 
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9. skew 


Purpose 


Skewness or third moment. 
Synopsis 
skewness=skew(x) 
D — 
skew calculates the skewness or third moment value. 
For vectors, skew(x) 1s the skewness value of the elements in vector x. 


For matrices, skew(x) is a row vector containing the skewness value of 


each column. 
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10. transfer 


Purpose 


Plot time series and frequency spectrum. 
Synopsis 
[Ya,y,fx]=transfer(yo,t,fo,nblock,sec) 
D ti 
transfer(yo,t,fo,nblock,sec) return three vectors, plot time senes and 


frequency spectrum. The return vectors are: 


Ya averaged fft vector of yo 
y calibrated velocity vector 
fx frequency vector 


The input parameters include: 


yo original amplitude 
t time series corresponds to yo 
fo operating frequency (Hz) 


nblock number of data blocks 
sec interested frequency sections 

The spectrum will indicate the peak position assign in sec. 
For example, input sec=[5 150], the program will search the peak point 


between 5 Hz and 150 Hz. 


a 


B. PROGRAM CODE 


1. kurt.m 


function kurtosis = kurt(y) 

% KURT Kurtosis. 

% For vectors, KURT(y) returns the kurtosis. 

% For matrices, KURT(y) 1s a row vector containing 
% the kurtosis of each column. See also COV. 


% Chao-Shih Liu 4-26-93 


[n,m]=size(y); 
if n==1 
vase. 
[n,m]=size(y),; 
end 
for1i=l:m 
meany (:,1)=mean(y(:,1)).*ones(n, 1); 
stdy(:,1)=std(y(:,1))*ones(n, 1); 
end 
s=((y-meany)./stdy ).%4; 
kurtosis = sum(s)/n, 
end 
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logbar.m 


function logbar( Vdb,L ineSty le,Line Width) 


% 

% logbar Plot 1/1 octave band bar plot using semilog in 

% x axis(freq.) and linear scale in Y axis(Vdb) 
% 

% logbar(Vdb) plots the semilog bar plot 

% for 1/1 octave band analysis. 

% 


% logbar(Vdb,LineStyle) plots the semilog bar plot 

% for 1/1 octave band analysis with various 

% types. See “plot” for more about LineStyle 

% 

% Line Width is a optional parameter to control the line 
% width of the bar. Default value is 0.5 

% 


% Wnitten By Liu,Chao-Shih 7/27/93 


if nargin<2 
lineStyle=’-’; 

end 

if nargin== 
linewidth=0.5, 

end 


% 

% Define the upper and lower limit of the pass bands (Hz) 
% 

flow=[11.2 22.4 45 90 180 355 710 1400 2800 5600], 
fhigh=[flow(2:10) 10000}, 


% 

% Determine the coordinates of the bars 
% 

xl=flow(:); 

xh=fhi gh(:); 

y=Vdb(:); 

[n,m]}=size(y); 

n3=n*3;, 

yy=zeros(n3+1,m); 

XX=yy; 

xx(1:3:n3,:)=x1(1 :n); 
xx(2:3:n3,:)=xh(1:n); 
xx(3:3:n3,:)=xh(1:n); 

xx=[xI(1);xx]; 

yy(1:3:n3,:)=y; 

yy(2:3:n3,:)=y; 

yy=[0;yy]; 
set(semilogx(xx,yy,LineStyle),’ linewidth’ Line Width) 
end 
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mstat.m 


function mstat(y,fs,nbar,mid,locat,manaly,mthre,z,mregen,pn) 
% 

% mstat(y,fs,nbar,mid,locat,manaly,mthre,z,mregen,pn) 
% will analysis the onginal data y and produce a 

% statistical report with a nbar pdf plot. 

% 

% mstat(y) will allow the user to perform the program in an 
% interactive way. 

% 

% y : amplitude (after calibration) 

% fs : sampling frequency (Hz) 

% nbar : number of bars in Prob. Density Distnbution Plot 
% mid : type of machinery 

% locat : pickup location 

% manaly: method of analysis 

% =] means All data points analysis 

% =2 means Overall peak points analysis 

% =3 means Cross peak points analysis 

% mthre : with threshold or not 

% =] without threshold 

% =2 with threshold 

% z: percentage of threshold 

% mregen: transform to Vdb or not 

% =] without transform 

% =2 transform to Vdb domain 

% pri: pmnt or not 

% =’y’ means to print out the result 


% 
% Liu,Chao-Shih 6/7/93 
% 


if nargin==1 

fs=input(‘Input the sampling frequency ? (Hz)’); 
nbar=input(‘Input the number of bar used in the PDF plot ?’); 
mid=input(‘Input the type of the machinery ?’,’s’); 
locat=input(‘Input the pickup location ?’,’s’); 
manaly=menu(‘Choose the method of analysis’... 
‘All data points analysis’,... 

‘Overall peak analysis’,... 

‘Cross peak analysis’); 

mthre=menu(‘Type of the threshold ?’,... 

‘Without threshold’... 

‘With n% threshold’); 

if mthre==2 

z=input(‘Input the % of the threshold’ );z=z/100; 
end 

mregen=menu(‘ Transform the data ?’,... 

‘Onginal data’,... 

‘Vdb Domain’), 
end 


% 
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% Method of analysis 
% 

if manaly==1 

ypry; . 
mt=’All data points’ 
elseif manaly== 
yp=npeak(y); 

mt=’ Overall peak’, 
else 

yp=npcly); 
mt=’Cross peak’; 
end 


% 

% Type of threshold 

% 

i mthre=— 

tht=’ without threshold’; 

else 

tht=[* with ©‘ num2str(z*100) ‘% threshold’); 
yp=nthre(yp,z), 

end 

pnumber=length(yp); 


% 

% Compute the first four mements for linear domain 
% 

yp=abs(y p), 

mepx=[‘m=’ num2str(mean(yp))]; 

stdpx=[‘s=’ num2str(std(yp))]; 
skewpx=[‘Skewness=’ num2str(skew(yp))]; 
kurtpx=[‘Kurtosis=’ num2str(kurt(yp))]; 


% 

% Compute the first four moments for Vdb Domain 
% 

if mregen==2 

yr=20* log! O(yp/le-8), 
meanyr=mean(yr);stdyr=std(yr); 
skewyr=skew(yr);kurtyr=kurt(yr); 
merpx=[‘m = ‘ num2str(mean(yr))]; 
stdrpx=[‘s = ‘ num2str(std(yr))]; 
skewrpx=[‘Skewness=’ num2str(skew(yr))]; 
kurtrpx=[‘Kurtosis=’ num2str(kurt(yr))]; 
end 


% 
% Plot the overall envelope PDF 
% 
clg, 

if mregen== 

subplot(312) 

[av,au }=hist(yp,nbar),bar(au,av); 
ylabel(‘Probability ’);xlabel(‘ Velocity (m/sec)’); 
title(“Onginal Data(x) Envelope’), 
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subplot(3 13);axis(‘off’); 

set(text(0.3,1,mepx), fontn’,’symbol’,’fonta’, italic’), 
set(text(0.3,0.8,stdpx), fontn’,’symbol’,’fonta’, italic’); 
set(text(0.3,0.6,skewpx), fontn’,’times’,’fonta’, italic’), 
set(text(0.3,0.4,kurtpx), fontn’,’times’,’fonta’, italic’); 
else 

subplot(323) 

[av,au]=hist(yp,nbar);bar(au,av), 

y label(‘Probability ’);xlabel(‘ Velocity (m/sec)’), 
title(‘(a) Linear Domain’); 

subplot(324) 

[av,auJ=hist(yr,nbar);bar(au,av); 

ylabel( ‘Probability ’);xlabel(‘ Vdb (0 Vdb=1e-8 m/sec)’); 
title(‘(b) Vdb Domain’), 

subplot(313);axis(‘off’), 
set(text(0.1,1,mepx),’fontn’,’symbol’,’fonta’, italic’), 
set(text(0. 1,0.8,stdpx),’fontn’,’symbol’,’fonta’, italic’); 
set(text(0. 1,0.6,skewpx), fontn’,’times’, ’fonta’,’ italic’), 
set(text(O. 1,0.4,kurtpx), fontn’,’times’,’fonta’, italic’); 
set(text(0.7,1,merpx), fontn’,’symbol’,’fonta’, italic’); 
set(text(0.7,0.8,stdrpx), fontn’,’symbol’,’fonta’, italic’); 
set(text(0.7,0.6,skewrpx), fontn’,’times’,’fonta’, italic’); 
set(text(0.7,0.4,kurtrpx), fontn’,’times’,’fonta’, italic’), 
end 


% 

% Statistical analysis results 

% 

subplot(311),axis(‘off’); 

tal=[‘Statistical Analysis Results for * mid]; 
set(text(0.5,1,tal),’fontsi’,14,’fontw’, bold’,’hon’,’c’), 
loc=[‘Pickup Location : * locat]; 

set(text(0.5,0.8,loc), fontsi’,14,’fontw’,’bold’,’hon’,’c’); 
ta2=[‘Data analysis by using : ‘ mt tht]; 
set(text(0.5,0.4,ta2),’hon’,’c’); 

block=[*‘Number of data points : “ num2str(pnumber)]; 
set(text(0.5,0.2,block),’hon’,’c’); 

frange=[‘The sampling frequency : ‘ num2str(fs) ‘ Hz’); 
set(text(0.5,0,frange), hon’,’c’); 


% 

% Print the graph or not 

% 

if nargin<4 

pn=input(‘Pnnt the plot ? (y/n)’,’s’); 
end 

if pn==’y’ 
set(1,’PaperPosition’,[1.25 2 6 8.5}) 
pnint -f] 

end 

end 
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npc.m 


function ycp=npc(y) 
% 


% ycp=npc(y ) will find the peak points between each two 
% zero-crossing. 

% 

% y : onginal data 

% ycp : cross peak data 

% 


% 
% Liu,Chao-Shih 4/26/93 
% 


y=y(); 
k=length(y); 
y=y-mean(y)*ones(k, 1); 
npc=0;]=1; 
ycp=zeros(k, 1); 
for n=2:k 
if sign(y(n)*y(n-1))==- 
h=n-1;npc=npct1; 
[yp,np]=max(abs(y (1:h))); 
np=np+l-1; 
yep(npc)=y (np); 
l=h+1; 
end 
end 
ycp=ycp(I:npc), 
en 


> 
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npeak.m 


function yp=mpeak(y ) 

% 

% yp=mpeak(y) will find the overall peak points. 
% 

% y : onginal data 

% yp : overall peak data 

% 


% 
% Liu,Chao-Shih 4/14/93 
% 


y=y(:); 

k=length(y),; 
y=y-mean(y)*ones(length(k), 1); 
yp=zeros(k, |); 

yp(1)=y(1); 

pce=l; 

for n=2:k-1 

nsl=n-1; 

nl=sign((y(n+1 )-y(n))*(y(n)-y(ns1))), 
if nl==-1 | ym)==y(ns1) 
pe=pc+l;yp(pe)=y(n), 

end 

end 

yp=yp(1:pe, 1); 

end 


82 


nthre.m 


function yt=nthre(y,z) 

% 

% yt=nthre(y,z) will remove the points with z% threshold 
% 

% y : onginal data 

% yt: data after threshold 

% z: percentage of the threshold 

% If z=1 then the program will not return the peak 

% peak data small than 1% of the maxinum peak value. 


% 

% Liu,Chao-Shih 4/14/93 
% 

y=sort(y(:)); 
ly=length(y), 
maxy=max(y); 
miny=z*0.01*maxy; 
for l=1:ly 

if y()>=miny 

po=l; 

yt=y(po:ly, 1); 

break 

end 

end 
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octave.m 


function oct=octave(y,11,12) 

% 

% OCTAVE Octave band. 

% 

% oct=OCTAVE(y ) retums the 1/1 octave band Vdb values from 
% band | to band 10. 

% 

% oct=OCTAVE(y,11,12) returns the 1/1 octave band Vdb values 
% from band 11 to band 12 and the others 

% bands will be assigned to 0. 

% 

% y(.,1): Vdb values 

% y(:,2) : frequencies 


% 
% Chao-Shih Liu 7/13/93 
% 


fc=[16 31.5 63 125 250 500 1000 2000 4000 8000]; 
flow=[11.2 22.4 45 90 180 355 710 1400 2800 5600}; 
fhigh=[22.4 45 90 180 355 710 1400 2800 5600 11200}; 
np=max(size(y )); 

oct=zeros(10,1); 


% 

% Default the 11 and 12 
% 

if nargin==1 
11=1;512=10: 

end 


% 
% if Vdb lower than 0 db, force the value=Odb 
% 

for k=l:np 

if y(k,1)<O 

y(k, 1)=0, 

end 

end 


% 

% convert to linear velocity value 
% 

vref=le-8, 
y(.,1)=10.%(v¢,1)/20)*vref; 
fmax=y(np,2),df=fmax/(np-1), 


% 

% Calculate the 1/1 octave band levels in linear domain 
% 

df2=df/2, 

fupper=y(:,2)+df2; % Upper limit of frequency 
flower=y(.,2)-df2; % Lower limit of frequency 
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sum=zeros(i2-11,1 ); 
for k=11:12; 
mark1=fupper-flow(k); 
mark2=flower-flow(k); 
mark3=fupper-fhigh(k), 
mark4=flower-fhigh(k); 
for j=1:np 
if mark4(j)>=0 % larger than the band(k) 
break 
else 
if mark] (})<=0 % smaller than the band(k) 
% do nothing 
else % cover the band(k) 
if mark3(j)>0 
if mark2(j)<O % band(k) smaller than df 
sum(k)=y(j,1)*2*(fhigh(k)-flow(k)); 
break; 
else % partial band in band(k+1) 
sum(k)=sum(k)+y(@, 1 )*2*abs(mark4(j)); 


end 
else 
if mark2(})<O % partial band in band(k-1) 
sum(k)=y(j,1)*2*mark1(j); 
else 
sum(k)=sum(k)+y Gj, 1)*2*df; 
end 
end 
end 
end 
end 


% Convert the octave band levels from linear scale to VdB levels 
oct(k)=20* log 1 0(sqrt(sum(k))/vref); 
end 
% Plot the results 
bar(oct) 
xlabel(‘ Number of Center Frequency ’) 
ylabel(‘ Vdb value (OVdb=1e-8 m/sec)’) 
axis({il-1 12+1 40 120]) 
for k=11 112 
octext=spnntf(°%S. If’ ,oct(k)); 
set(text(k ,oct(k)+2,octext), ’fontw’, bold’, ’horiz’, center’) 
end 
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psd.m 


function Vdb=psd(y,nblock) 

% Vdb=PSD(y,nblock) 

% 

% This function will compute the average Vdb value by 
% using hanning to nblock data set. 

% 

% y : original data matnix (time,velocity) 

% nblock : number of data blocks 

% Ya: average fft Vdb vector of y(OVdb=1le-8 m/sec) 
% fx : frequency vector 

% 

% Liu,Chao-Shih 12:11PM 7/14/93 

% 


% standardize the data matnx 
n=size(y,2); 

if n~=2 

ay 

end 


t=yC, Dsy=yC2), 
dt=t(2)-td1 ); 

fs=1 /dt, 
vref=1e-8. 


% 
% Compute the fft data 
% 
n=length(y )/nblock;n2=n/2; 
nl=1-; 
w=hanning(n); 
for nb=] :nblock 
nh=nb*n; 
Y(1:n,nb)=fft(w. *y (nl:nh,1))*dt; 
nl=nh+1 ; 
end 
if nblock~=1 
Y=(sum(abs(Y )’)/nblock)’; 
else 
Y=abs(Y), 
end 
Ya=20*log10(2* Y/vref); 
f=fs/2*(0:n2)/n2; 
Ya(n2+2:n)=[]; 
fx=f(2:n2+1)’; 
df=fx(2)-fx(1); 
fx=[0;fx]; 
Vdb=[ Ya, fx]; 
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skew.m 


function skewness = skew(y) 

%SKEW Skewness. For vectors, SKEW(y) returns the skewness. 
% For matrices, STD(y) is a row vector containing 

% the skewness of each column. See also COV. 


% Chao-Shih Liu 4-26-93 


[n,m]=size(y), 
if n== 

yy; 
[n,m]=size(y ); 
end 
for i=1l:m 


meany(:,i)=mean(y(:,1)).*ones(n, 1); 
stdy(:,1)=std(y(:,1))*ones(n, 1); 
end 
s=((y-meany)./stdy ).*3; 
skewness = sum(s)/n; 
end 


87 


10. transfer.m 


function [ Ya,y,fx }=transfer(yo,t,fo,nblock,sec) 
% [Ya,y,fx]=transfer(yo,t,fo,nblock,sec) 

% 
% This function will plot time & frequency domain 
% 

% yo : onginal velocity vector (before calibration) 
% t : time vector 

% fs : sample frequency 

% fo : operation frequency = RPM/60 

% nblock : number of data blocks 

% Ya: average fft vector of y 

% fx : frequency vector 

% y : onginal velocity vector (after calibration) 

% 

% Liu,Chao-Shih 1:58 7/10/93 

% 


dt=t(2)-t(1); 

fs=1/dt; 

itext=input(‘Input the title of this record =>’,’s’); 

locat=input(‘Input the location of data collect’,’s’); 

y=yo, 

scale=menu(‘Choose the method of X scale’,’Frequency (Hz)’, Order’); 
vref=1e-8; 


% 

% Standarized the size of y 
% 

[xd,yd]=size(y); 

if yd~—=1 

ay et 

end 

ny=length(y),; 


% 

% Compute the fft data 

% 

n=ny/nblock;n2=n/2; 

nl=1; 

win=menu(‘Type of Window System’,’ Without window system’,... 
‘Hanning window’,’ Hamming window’), 
if win== 

w=ones(n, |); 

elseif win==2 

w=hanning(n),; 

else 

w=hamming(n), 

end 

for nb=1:nblock 

nh=nb*n, 

Y(1:n,nb)=fft(w. *y(nl:nh,1))*dt; 
nl=nh+1-; 
end 
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if nblock~=1 

Y=sum(abs(Y)’)/nblock; 

else 

Y=abs(Y); 

end 

Ya=20* log] 0(2* Y/vref); 

f=fs/2*(0:n2)/n2; 

Ya(n2+2:n)=[]; 

if scale==1 

fx=f(2:n2+1); 

df=fx(2)-fx(1); 

xt=’Frequency (Hz)’;unit=’ Hz ‘; 

else 

fx=f(2:n2+1 )/fo; 

df=fx(2)-fx(1); 
low=ceil(0.75/df);high=ceil(1.2S/df); 

[Ymax, Ypos]=max(Ya(low:high));fscale=fx(low+ Ypos-1),; 
fx=fx/fscale;, fo=fo*fscale; 

xt=[‘Frequency (Normalized Orders)’];unit=’ Order ‘; 
end 
df=fx(2)-fx(1); 

sec=ceil(sec/df),; 

if sec(1)==0 

sec(1)=1; 

end 

for mk=]1 :length(sec)-1 
[mav(mk),pos(mk)]=max(Ya(sec(mk):sec(mk+1 ))); 
pos(mk)=pos(mk)+sec(mk), 

end 


% 
% Plot time domain and frequency domain 
% 
figure(1);clf;subplot(2 1 2), 
tf=1024; 
if ny<1024 
=ny; 

end 
plot(td1:tf,1),y (1:tf, 1),’r’);axis((t(1,1) t¢tf,1) min(y) max(y)]); 
title(‘ Time Domain’);xlabel(‘Time (sec)’);ylabel(‘ Velocity(m/sec)’); grid 
% 
% General information 
% 
subplot(211);axis(‘off’) 
mfo=(‘The operation condition = ‘ num2str(fo*60)... 

‘RPM ©‘ num2str(fo) ‘ Hz)’]; 
mloc=[‘Data pickup location = ‘ locat]; 
mfs=[‘Data sampling frequency = ‘ num2str(fs) ‘ Hz’); 
mdp=[‘Total number of sampling data points = ‘ num2str(length(y))]; 
set(text(0.5,],1text),’fontsi’,15,’fontw’,’bold’, horiz’,’center’) 
text(0.05,0.9,mloc);text(0.05,0.82,mfo); 
text(0.05,0.74,mfs);text(0.05,0.66,mdp); 


text(0,0.5,’* The time domain’), 
mt=(‘Mean value = ‘ num2str(mean(y)) ‘ m/sec’); 
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text(0.05,0.4,mt); 

st=[‘Standard deviation = ‘ num2str(std(y)) ‘ m/sec’); 
text(0.05,0.32,st); 

text(0.05,0.24,[‘The plot only shown * num2str(tf) * data points’); 
pause 


figure(2);,clf; subplot(2 1 2), 

plot(fx, Ya(2:n2+1),’r’); 

upper=(ceil((max(mav))/10)+1)* 10; 

axis({O fx(n2) 20 upper]), 

title(‘Frequency Domain’);xlabel(xt);y label(‘ Velocity in Vdb’);gnd 
for lop=1:length(sec)-1 

text(fx(pos(op)-1),mav(lop),num 2str(lop)); 

end 


% 

% General information 

% 

subplot(211);axis(‘off’) 

mfo=[‘The operation condition = ‘ num2str(fo*60)... 

‘RPM (= ‘ num 2str(fo) ‘ Hz)’); 

mloc=[‘Data pickup location = * locat]; 

mfs=[‘Data sampling frequency = ‘ num2str(fs) ‘ Hz’); 
mbk=[‘Number of data block averages = * num2str(nblock)]; 
mdp=[‘Block size = ‘ num2str(length(y )/nblock) ‘ points/block’}; 
set(text(0.5,1,itext), fontsi’,15,’fontw’,’bold’,’honz’,’center’),; 
text(0.05,0.9,mloc);text(0.05,0.82,m fo); 
text(0.05,0.74,mfs);text(0.05,0.66,mbk);text(0.05,0.58,mdp), 
text(0,0.5,’* The frequency domain’), 

for lop=1 :length(sec)-1 

as=[num2str(lop),” = (‘,num2str(fx(pos(lop)-1)).... 

unit,’, ‘,.num2str(mav(lop)),’ db)’); 

if lop>4 

xpos=0.55;ypos=0.8-0.08* lop; 

else 

xpos=0.05;ypos=0.48-0.08* lop; 

end 

text(xpos,ypos,as); 

end 

morder=[‘l order = ‘ num2str(fo) * Hz’}; 

text(0.05,0.08 morder),; 

mref=[‘The reference 0 Vdb = ‘ num2str(vref) ‘ m/sec’); 
text(0.05,0,mref); 

% 

% zooming up the plot ? 

% 

zoom= 1. 

while zoom~=3 

subplot(212) 

zoom=menu(‘zoom up ?’,’0-10 order’... 

‘Use mouse to select the range’,’No’); 

if Zoom== 

axis({O 10 20 upper]); 

elseif zoom== 

axis((O fx(n2) 20 upper]); 
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[xs,ys]=ginput(2); 
xrange=sort(xs); 

if xrange(1)<O 

xrange(1)=0; 

end 

axis({xrange(]) xrange(2) 20 upper]); 

else 

break 

end 

end 

% 

% Pnint out the plot 

% 

uiprint=[ ‘set(1,”’PaperPosition’’,[0.25 1 8 9.25]);’ ‘print -f1,’ ... 
‘set(2,”’PaperPosition”’,[0.25 1 8 9.25]);” ‘print -f2;’ ‘end’); 

uicontrol(‘style’, ’push’,’units’,’normal’,’pos’,[0.9 0.9 0.08 0.06], ... 
‘string’, ’print’,’call’ ,uipnint), 
uicontrol(‘style’,’push’,’units’,’normal’,’pos’,[0.9 0.82 0.08 0.06], ... 
‘string’, ’exit’,’call’,’end’); 
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APPENDIX B. EXAMPLES OF USING MATLAB PROGRAM IN 
PC486 


MATLAB allows the user to wmite a script file. A script file is an external file that 
contains a sequence of MATLAB statements. The following example will show the script 


files used to perform the analysis. 


A. CROSS PEAK ANALYSIS EXAMPLE 


In this example, a CG-59 fire pump measured at MB(FE) data set is used in the 
analysis. The file name of the data file is f/p/.asc. The MATLAB command is anchored 
by >>. 

>> load flp1.asc 

>> y=flp1(:,2); 

>> t=flpl(:,1); 

>> mstat(y) 


The user can then interact with the program to select the method of analysis. Instead 
of this interactive process, the user can use the following to execute the program in an 
automatic way. 

>> load f1pl.asc 

>> y=tlpl(:,2); 

>> t=flpl(.,1); 

>> ypc=npc(y); 

>> yt=nthre(ype, 1); 

>> mid=’CG-59 Fire Pump’; 

>> locat=’MB(FE)’; 
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>> mstat(yt, 10000, 10,2,mid,locat,3,2,1,2,’y’) 


B. 1/1 OCTAVE BAND ANALYSIS EXAMPLE 

Here, proceeding as in the above example, a CG-59 fire pump measured at MB(FE) 
data set is used to perform the 1/1 octave band analysis. The file name of the data file is 
fipl.asc. The MATLAB command 1s anchored by >>. 


>> load flp1].asc 


>> yh=f1p1(:,2); ‘Yo yh is a time series sampled at 10 kHz. 
>> yl=yh(1:5:length(yh)); Yo yl is a time series sampled at 2 kHz. 
>> t=flp1(,1); 


>> VdBh=psd(yh,20); 

>> VdBl=psd(yl,20); 

>> octh=octave(VdBh,7,9); 
>> octl=octave(Vdbl, 1,6); 
>> octband=[octh octl]; 


>> logbar(octband); 
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APPENDIX C. FIGURES OF STATISTICAL ANALYSIS RESULTS 


This section contains figures of the statistical analysis results with and without 


considering damaged data sets. 


94 


Statistical Analysis Results 
Type | Fire Pump 
Pickup Location : MB(FE) 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 1.965e+005 
The sampling frequency : 1e+004 Hz 


(a) Linear Domain (b) VGB Domain 
0.6 0.2 
204 2 
2 2 
2 @ 0.1 
© 2 
oO. 0.2 a: 
0 0 
0 2 4 60 80 100 120 
Velocity (m/sec) x 10° VdB (OVdB=18-8 nvsec) 
u=0.0004954 pL = 89.36 
o=0.0005325 o= 9.203 


Skewness=1.85 Skewness=-0.04575 


Kurtosis=6.201 Kurtosis=2.378 


Figure C.1 The Statistical Analysis Results for Type I Fire Pumps Measured at 
MB(FE) (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type | Fire Pump 
Pickup Location : MB(CE) 


Data analysis by using : Cross peak with 1% threshold 


Number of data points : 2.173e+005 


The sampling frequency : 1e+004 Hz 


(a) Linear Domain 
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Skewness=2.052 
Kurtosis=7.195 


(b) VadB Domain 
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oO 
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pt = 90.02 
o=9.12] 
Skewness=0.04893 
Kurtosis=2.444 


Figure C.2 The Statistical Analysis Results for Type I Fire Pumps Measured at 
MB(CE) (a) Linear Velocity Domain (b) VdB Domain 


Statistical Analysis Results 
Type | Fire Pump 
Pickup Location : PB(CE) 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 3.129e+005 
The sampling frequency : 1e+004 Hz 





(a) Linear Domain (b) VdB Domain 
0.6 0.2 

20.4 = 

2 5 

& 80.1 
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Velocity (m/sec) VdB (OVdB=1e-8 m/sec) 
u=0.002828 pL = 104.6 
o=0.00268 Oa 7) 
Skewness=1.517 Skewness=-0.4436 
Kurtosis=5.638 Kurtosis=2.476 


Figure C.3 The Statistical Analysis Results for Type I Fire Pumps Measured at 
PB(CE) (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type | Fire Pump 
Pickup Location : PB(FE) 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 2.614e+005 
The sampling frequency : 1e+004 Hz 





(a) Linear Domain (b) VdB Domain 
0.6 
20.4 = 
a TaN 
wo oO 
Q ai 
o o 
o 0.2 as 
0 0 —— ) 

0 0.005 0.01 0.015 50 100 150 
Velocity (m/sec) VdB (OVdB=16-8 nvsec) 
u=0.002717 pt = 103.6 
o=0.00268 1 C= 10/45 
Skewness= 1.253 Skewness=-0.3354 
Kurtosis=3.959 Kurtosis=2.232 


Figure C.4 The Statistical Analysis Results for Type I Fire Pumps Measured at 
PB(FE) (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type | Fire Pump 
Pickup Location : PB(FE/A) 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 3.105e+005 
The sampling frequency : 1e+004 Hz 


(a) Linear Domain 





(b) VdB Domain 
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Figure C.5 The Statistical Analysis Results for Type I Fire Pumps Measured at 
PB(FE/A) (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type Il Fire Pump 
Pickup Location : MB(FE) 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 4.733e+004 
The sampling frequency : 1e+004 Hz 


(a) Linear Domain (b) VdB Domain 
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Figure C.6 The Statistical Analysis Results for Type Il Fire Pumps Measured at 
MB(FE) (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type Il Fire Pump 
Pickup Location : MB(CE) 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 4.835e+004 
The sampling frequency : 1e+004 Hz 





(a) Linear Domain (b) VdB Domain 
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Figure C.7 The Statistical Analysis Results for Type II Fire Pumps Measured at 
MB(CE) (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type Il Fire Pump 
Pickup Location : MB(CE/A) 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 3.089e+004 
The sampling frequency : 1e+004 Hz 


(a) Linear Domain (6) VdB Domain 
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Figure C.8 The Statistical Analysis Results for Type II Fire Pumps Measured at 
MB(CE/A) (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type Il Fire Pump 
Pickup Location : PB(CE) 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 3.241e+004 
The sampling frequency : 1e+004 Hz 


(a) Linear Domain (b) VdB Domain 
0.6 0.2 
210.4 = 
a ‘o 
xs aS 0.1 
2 o 
ao 0.2 o 
0 0 
0 0.005 0.01 60 80 400 120 
Velocity (m/sec) VdB (OVdB=1e-8 m/sec) 
u=0.002484 u= 102.6 
o=0.002344 o= 10.84 
Skewness=0. 9047 Skewness=-0.3675 
Kurtosis=2.675 Kurtosis=2.038 


Figure C.9 The Statistical Analysis Results for Type I Fire Pumps Measured at 
PB(CE) (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type II Fire Pump 
Pickup Location : PB(FE) 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 4.056e+004 
The sampling frequency : 1e+004 Hz 


(a) Linear Domain (b) VdB Domain 
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Figure C.10 The Statistical Analysis Results for Type If Fire Pumps Measured at 
PB(FE) (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type Il Fire Pump 
Pickup Location : UMB 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 2.1e+005 
The sampling frequency : 1e+004 Hz 


(a) Linear Domain (b) VdB Domain 
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Figure C.11 The Statistical Analysis Results for Type III Fire Pumps Measured at 
UMB (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type Ill Fire Pump 
Pickup Location : LMB 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 2.176e+005 
The sampling frequency : 1e+004 Hz 
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Figure C.12 The Statistical Analysis Results for Type III Fire Pumps Measured at 
LMB (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type Ill Fire Pump 
Pickup Location : LMB(A) 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 2.845e+005 
The sampling frequency : 1e+004 Hz 


(a) Linear Domain (b) VdB Domain 
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Figure C.13 The Statistical Analysis Results for Type III Fire Pumps Measured at 
LMB(A) (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type Ill Fire Pump 
Pickup Location : UPB 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 1.367e+005 
The sampling frequency : 1e+004 Hz 


(a) Linear Domain (b) VdB Domain 
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Figure C.14 The Statistical Analysis Results for Type III Fire Pumps Measured at 
UPB (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type Ill Fire Pump 
Pickup Location : LPB 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 1.99e+005 
The sampling frequency : 16¢+004 Hz 
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Figure C.15 The Statistical Analysis Results for Type III Fire Pumps Measured at 
LPB (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type | Fire Pump 
Pickup Location : PB(CE) w/o 89#6,90#3,90#6 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 2.704e+005 
The sampling frequency : 1e+004 Hz 


(a) Linear Domain (b) VdB Domain 
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Figure C.16 The Statistical Analysis Results (Without Damaged) for Type I Fire 
Pumps Measured at PB(CE) (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type | Fire Pump 
Pickup Location : PB(FE) w/o 90#3 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 2.489e+005 
The sampling frequency : 1e+004 Hz 
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Figure C.17 The Statistical Analysis Results (Without Damaged) for Type I Fire 
Pumps Measured at PB(FE) (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type Il Fire Pump 
Pickup Location : MB(CE) w/o #1 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 3.238e+004 
The sampling frequency : 1e+004 Hz 
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Figure C.18 The Statistical Analysis Results (Without Damaged) for Type II Fire 
Pumps Measured at MB(CE) (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type Ill Fire Pump 
Pickup Location : UMB w/o 11/24/#1 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 2.05e+005 
The sampling frequency : 1e+004 Hz 
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Figure C.19 The Statistical Analysis Results (Without Damaged) for Type II Fire 
Pumps Measured at UMB (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type Ill Fire Pump 
Pickup Location : LMB w/o 11/24/#1,3/10/#2 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 1.873e+005 
The sampling frequency : 1e+004 Hz 


(a) Linear Domain (b) VdB Domain 
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Figure C.20 The Statistical Analysis Results (Without Damaged) for Type III Fire 
Pumps Measured at LMB (a) Linear Velocity Domain (b) VdB Domain 
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Statistical Analysis Results 
Type Ill Fire Pump 
Pickup Location : UPB w/o 11/24/#1 
Data analysis by using : Cross peak with 1% threshold 
Number of data points : 1.245e+005 
The sampling frequency : 1e+004 Hz 
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Figure C.21 The Statistical Analysis Results (Without Damaged) for Type III Fire 
Pumps Measured at UPB (a) Linear Velocity Domain (b) VdB Domain 
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